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We consider the influence of magnetic excitations on the local density of states in the d-wave 
superconductor. The magnetic susceptibility is calculated within the renormalized t — t — J model 
and its influence on the quasip article self-energy i s considered using a minimal model originally 
proposed by Polkovnikov et al. [ :ond-mat / 0203176{ . We find the local density of states possess 
periodic components both along (n, 0) and (tt, it) directions with the associated wavevectors changing 
in magnitude as the quasiparticle energy is varied. Comparison with the STM experiment reveals 
that the calculated LDOS modulation is inconsistent with the measured data. 



Introduction: Observation of discernible "checker- 
board" patterns in the local density of states 
(LDOS) has been reported for a superconducting 
Bi2Sr2CaCu208+5 (BSCCO) compound in the mixed 
state, with the spatial periodicity close to four lattice 
constants Q. Howald et al. reported the same peri- 
odicity exists for the LDOS in the same compound but 
without the magnetic field Hoffmann et al. in turn 
argued that this periodicity is in fact a function of energy, 
i.e. as the probing bias of the scanning tunnelling micro- 
scope (STM) is varied the spatial modulation periodicity 
also changes systematically ||. 

A number of theories have emerged in response to such 
striking set of observations . Some of these theories 
emphasize the possible emergence of a new order param- 
eter near the vortex core Q but, in view of the latest ex- 
perimental finding || may have difficulty to explain how 
an induced order parameter can exhibit a length scale 
that depends smoothly on energy. An alternative expla- 
nation advanced by Wang and Lee (WL) ascribes the pe- 
riodicity to the quasiparticle scattering due to localized 
impurities || . In this approach scattering processes con- 
nect parts of the underlying Fermi contours at different 
wavevectors as the quasiparticle energy is varied. 

Quite distinct from the proposal of WL is the theory of 
Polkovnikov, Vojta, and Sachdev(PVS) © who showed 
in a simple model that a dynamic spin fluctuation mode 
with a distinct ordering wavevector Q may lead to ob- 
servable periodic modulations in the local STM spectra 
at twice Q. They also invoke the local impurity pinning, 
which is necessary to break the translational symmetry 
and give rise to spatially varying LDOS. In principle the 
mechanism proposed by PVS does not depend on the 
presence of an external magnetic field. 

Dynamical spin fluctuations are observable in a num- 
ber of families of cuprates. It has been shown that 
magnetic absorption peaks occur at incommensurate 
wavevectors (tt, tt ± 8), and (tt ± 8, 7r) with 8 varying with 
the absorption energy [Q. Brinckmann and Lee (BL) 
Q argued that the underlying Fermi surface topology in 
the superconducting state is responsible for such incom- 
mensurate response. Within the Gutzwiller-renormalized 



t—t' — J model they indeed found the absorption peaks at 
(tt ± 8, tt) , and (tt, tt ± 8 ) while at even lower energies, the 
incommensurate peaks appear in the diagonal direction, 
(tt±8',tt±8'). 

Coupled with the idea of PVS, it is worthwhile to ask 
whether such incommensurate magnetic spectrum will 
also modulate the LDOS at wavevectors that vary with 
the energy. In this paper, we examine this issue using 
the magnetic spectrum calculated from the renormalized 
t — t' — J model. Assuming a simple interaction scheme 
between the spin of the quasiparticles and the mag- 
netic excitation we evaluate the quasiparticle self-energy 
and thence the spatially varying LDOS. Ultimately our 
calculation is aimed at a comparison with the latest ex- 
perimental data || , and in the process hope to figure out 
to what extent, if at all, the magnetic scenario for the 
observed periodicity is viable. 

Methods: We adopt the renormalized t — t' — J model 
H as the starting point: 
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In the above, t eff = tx + (J/^MjMXff = 
t'x,Aij = (e a pfiafj/3), and x is the average doping. 
J e ff is the residual spin-spin interaction between quasi- 
particles. 

From the above mean-field Hamiltonian one calculates 
the zero-temperature magnetic susceptibility xo(q ; * l/ )i 
and through the residual spin-spin interaction J e ffSi -Sj, 
the renormalized spin susceptibility 



x(q ) " / ) = 



xo(q,w) 



i + .-/(q)xo(q, iv) 



(2) 



for J(q) = J e f f (cos q x +cosq y ). In PVS's treatment the 
magnetic susceptibility is approximated by a phenomeno- 
logicalform, Xpvs(q, iv)~ x oc v 2 +A^+c 2 (q— Q) 2 , where 
A s is the spin gap, and Q is the ordering wavevector. 
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On the other hand, the t — t' — J model result for the 
renormalized magnetic susceptibility is known to exhibit 
a continuous set of ordering wavevectors with respect to 
energy S . We show a horizontal scan of the imaginary 
part x"(q.) v ) with q y = ir for a series of energies u> in Fig. 
|l|. The incommensurate peaks near q x — ir are obvious, 
and both these peaks and the less pronounced peaks near 
q x = move toward q x = ir when the energy is raised. In 
particular the peaks at (ir ± 5(E), ir) may give rise to a 
modulation in the charge sector with twice the wavevec- 
tor @§, i.e. at (2tt ± 25(E), 2ir) ee (±25(E),0). With 
5(E) being a decreasing function of energy, one naively 
expects LDOS modulation wavevectors to decrease at a 
higher bias as well, in qualitative agreement with the 
experiment Below, we examine through explicit cal- 
culation whether this naive argument will hold. 
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FIG. 1. Imaginary part of the magnetic susceptibility, 
x"(q, v) for q y — ir, and < q x < 2ir calculated from the 
renormalized t — t' — J model with J e // = 0.34J. Different 
plots correspond to v = 0.12J, 0.16J, 0.20J, and 0.24J. 

The treatment can proceed by positing the effective 
action C = Cm + C qp + C qv -u consisting respectively 
of the magnetic and quasiparticle sectors and their cou- 
pling. For the magnetic part ||, 

C M [®,®*] = \ J d 3 xd 3 y^(x) X - 1 (x-y)^(y) 

-IJdTlC&iTirt+WiruT)]. (3) 

Space-time coordinates are denoted by x — (r, r) and 
y = (r',r'). The last term is the pinning potential, lo- 
cated around r,;, for the spin density wave(SDW). Due to 
the pinning term the propagator ($*(r, iuj)&(r' , iuj)) = 
\(y,y' ,iui) is no longer translationally invariant. 

The quasiparticle coupling to the fluctuating SDW is 
modelled by § 

£ qp .M = \gS(vr) ■ [$(rr) + (rr)] . (4) 

Through this term the quasiparticle propagator is influ- 
enced by the collective mode, and exhibits translational 
symmetry breaking. The perturbative treatment of the 
effective interaction is straightforward, with the one-loop 
result for the change in Green's function 5Q pi| . 



5G(r,r',iu)= ^ 9 2 { iv )x( v i- v iiiv)T(iv)x(Yi-r 2 ,w) 

x G(r—r 1 ,iw)G(ri—r 2 ,iuj — w)G(r2 — r / ,iuj), (5) 

where T(iv) = x(r = 0, if). Frequency dependence in 
the coupling g(iv) is introduced for generality. In ar- 
riving at this expression we have only kept those terms 
that break the translational symmetry and ignored the 
space-independent self-energy correction. After the Wick 
rotation, iuj — > u> + i5, the spatially varying LDOS is 
N(r,uj) = (l/7r)Im<%?(r,r,w + iS). 

Normally the summation Y]^, is carried out analyt- 
ically, and the resulting expression is evaluated by nu- 
merical means after the Wick rotation ioj — > lo + i5. In 
our case the renormalized susceptibility x(Q: w) does not 
possess a simple Lehmann expansion, which prevents the 
frequency summation from being carried out. Instead 
we resort to the time-ordered Green's function method 
(as opposed to the Matsubara Green's function method) 
which does not require the Wick rotation. At zero tem- 
perature LDOS can be obtained from the time-ordered 
Green's function by (T =time-ordered) 

N(r, lu) = (l/7r)sgn(w)Im<5S T (r, r, u). (6) 

Perturbation theory proceeds in entirely analogous man- 
ner, and we only quote the final result for the Green's 
function correction. 

5G t (t,t',uj) = -i ^ / 9 2 {v)XT(-ri-Y l ,v)T(v)xT(-r l -v 2 ,v) 

x Gr(r— y\,uS)Gt (yi—t^w-v) G r (r2-r',w). (7) 

The time-ordered susceptibility xt and Green's function 
Gt are obtained from their Matsubara counterparts by 
iuj — > lu + isgn(w). 

Results: We evaluate Eq. ((?]) numerically on a large 
lattice of dimension [201x201] using self-consistently de- 
termined parameters A max = 2\(e a j3f ia fjp)\ = 0.48 J. 
The impurity site is located at the center, = (101, 101). 
We use J e ff = 0.34 J, and the broadening factor 
6 = 0.01J in evaluating x(r, v) and G(y,lo). The v- 
integration is replaced with a discrete sum, J dvf(v) — > 
& v Hs f( s ^ v )-> where sis an integer, -s max < s < s rnax . 
The following results are taken with Ai/ = 0.04J, and 
Smax = 14 implying that virtual processes with the off- 
shell energy of up to 0.56 J are taken into account in the 
self-energy evaluation. The commensurate response at 
q = (ir, 7r) in the susceptibility occurs at v res w 0.32 J 
for the present model, while at energies above and below 
v res , the resonance occurs away from (ir, it) fa. Our self- 
energy energy calculation therefore includes effects from 
both above and below the resonance energy. We have 
confirmed on a smaller lattice that using a finer scale 
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Av = 0.02 J does not change the final result qualitatively. 
The on-site Green's function 6Qt(t, r, oS) is evaluated for 
a subset of the lattice, typically for r lying in [101x101] 
about the impurity center, for \u>\ up to 0.32J. The real- 
space LDOS thus obtained is then Fourier-transformed 



according to 2V(q, uj) = \ J2 r 



c N(r,uj)\. 



FIG. 2. Fourier-transformed LDOS iV(q, uj) for the q-range 
[— 7r, 7r] x [n, 7r] and for —0.32 J < oj < 0.32 J. Center of each 
rectangle corresponds to q = (0,0). The intensity scheme 
is with respect to each figure, and should not be compared 
between different plots. (Please refer to attached figure 
Idoskl00b-v2.jpg) 

Plots in Fig. show iV(q, uj) calculated with the cou- 
pling constant g 2 ^) =const. Apart from the q = 
peak, we find some prominent peaks along both diago- 
nal (meaning (±7r,±7r)) and horizontal (meaning (±7r,0) 
and (0, ±7r)) directions. The peaks are located near the 
ends of the Brillouin zone at a small bias, and move to- 
ward q = as the bias energy is increased. A hori- 
zontal(diagonal) scan of iV(q, uj) allows the identification 
of wavevectors q_h(u->) (qd(uj)) for which iV(q, oj) is a lo- 
cal maximum. The values of qd{<-j) = |qd(^)|/v / 2 and 
qh{u) = |q/i(^)| are plotted in Fig. f|(a). Clearly both 
qd{u) and qh(u>) are decreasing functions of energy. The 
detailed shape of iV(q, uj) however depends somewhat on 
whether V is negative (electrons tunneling out) or posi- 
tive (electrons tunneling in). 

As shown in Fig. |l| the magnetic absorption peaks 
are most pronounced around (-7T, 7t). Note however that 
the Green's function shown in Eqs. (g) and (0) are a 
product of several complicated functions, and the mag- 
netic absorption spectrum alone does not a priori char- 
acterize the LDOS. In this regard we can loosely classify 
the structures in x"(q, v) as those near (it, 0), (0, it), and 
those near (ir, n). In the result shown in Fig. || con- 
tributions from these two regions will undoubtedly be 
mixed. Instead, if we truncate the magnetic spectrum 
to be within a certain vicinity of (tt, it) , the quasiparticle 
Green's function will also be influenced by the magnetic 
excitations at these wavevectors, and none from around 
("7T, 0), (0, 7r). This is achieved by using the reduced sus- 
ceptibility Xred(q, v) = x(q, v) for -7r/4 < q x , q y < tt/4, 
Xred(q, v) — otherwise. Its Fourier transform Xred(r, v) 
is used in the evaluation of the new LDOS, N re d(r,uj) 
and N re d((l, lu) , shown in Fig. In this calculation 
we use the cutoff g 2 (v) — g 2 {0) exp[— 0.1(^/Ai/) 2 ]. For 
a less sharp cutoff such as g 2 {y) = g 2 (0) or g 2 (v) = 
g 2 (0) exp[— 0.1|z^/AH], we nn d no other distinguishable 
structure in A^ re( ;(q, uj) except a broad hump centered at 
q = 0. 



FIG. 3. LDOS iV re( j(q, uj) calculated with the reduced mag- 
netic susceptibility \red- Energy and momentum ranges are 
the same as in Fig. ^| (Please refer to attached figure 
ldosk200e.jpg) 

There is very little overlap between the two figures 
§ and [| as also reflected in the qd{u) and qh{oj) val- 
ues obtained respectively from each figure. Figure ||(b) 
shows the peak intensity position of N re d{q_,oj). Essen- 
tially all qd(oj), qh{w) are greater than 7r/2 in the iV(q, uj) 
while in the reduced LDOS, all qd(oj), and <j7i(w) are less 
than 7r/2. Furthermore, qd(^) and qh{u) in Fig. [|(b) 
are increasing functions of energy, as opposed to Fig. 
|](a) where they are decreasing. In both plots, qd{^) 
and qh{u) show the same dependence on energy, being 
both increasing or both decreasing functions. As previ- 
ously mentioned a standard Ginzburg-Landau argument 
about the coupling between SDW and charge-density- 
wave (CDW) H implies that an ordering tendency occur- 
ring at (vr± 5, tt± 5) for the SDW will lead to CDW mod- 
ulations at (±28,±2S). Naively, since S decreases with 
energy, one expects qh(u) in Fig. ||(b) to be a decreasing 
function of uj, which is contradicted by our calculation. 
In obtaining Fig. [|we (1) truncated the magnetic spec- 
trum within a small momentum window around (-7T, tt) 
and (2) modified the coupling g 2 {v). With this manip- 
ulation we were able to bring out features in the LDOS 
which were "hidden" in the calculation that led to Fig. 

Since the coupling constant g must in reality be a 
function of both momentum and energy, one cannot be 
too certain a priori which of the features found in Figs. 
| and | are more readily observable in an experiment. 





FIG. 4. Energy dependence of qd{uj) and qh(uj) (defined 
in the text) obtained from JV(q, uj) (Figure (a)) and from 
N re d(<\, uj) (figure (b)). In each figure, four data sets are 
shown for q d (uJ < 0) (nd), q h (uj < 0) (nh), q d (uJ > 0)(pd), 
and qh(io > 0)(ph) against \oj\. Note that all the q values 
decrease with energy \uj\ in (a), and increase in (b). 

Discussion: We now come to comparison of our re- 
sult with the experimental data. First, it has been found 
in Ref. |^] that while 5d(w) increases, qh{u) decreases 
at a larger bias |u;|. This is in contrast to our calcu- 
lation where both qdi^j) and <?d(w) behave analogously 
with respect to energy. Secondly, experimentally mea- 
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sured kh(u>) values are found to lie below it/ 2 for a wide 
range of energies and hole doping, which may be consis- 
tent with Figure ||(b), but in this case the energy depen- 
dence is opposite to our calculation. In Fig. [|(a), <7d( w ) 
and qh(u) values are rather high which, translated into 
real space, implies LDOS modulations on the length of 
2 — 4 lattice spacings. While such magnetic-fluctuation- 
induced LDOS variations are not forbidden and may very 
well be observable in the future, we must conclude that 
the LDOS modulation based on the t — t' — J model of 
the magnetic excitation is inconsistent with the currently 
known experiment S. 

It is possible that the renormalized t — t' — J model does 
not after all capture the quasiparticle band structure of 
the BSCCO and that in another model, one indeed finds 
the energy dependence of qd and qh consistent with the 
experiment within the magnetic fluctuation scenario. For 
example, in WL Q , the model adopted is the phenomeno- 
logical tight-binding model originally proposed by Nor- 
man jyj. 

On the theoretical side, one should be careful that the 
naive Ginzburg-Landau argument for SDW-CDW cou- 
pling may in some cases lead to predictions which are in- 
consistent with a full many-body calculation such as this 
one. Although the patterns in iV(q, us) are ultimately 
due to the underlying magnetic fluctuation, the reason 
for the particular energy dependence of the modulation 
period found in Figs. ^ and [| remains unclear. Finally 
an entirely different mechanism such as WL's are not in- 
compatible with the present model, and may well simul- 
taneously lead to observable effects in a given system. 
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